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ABSTRACT 

Introduction: Migratory birds act as hosts for a number of zoonotic viruses, and have the 
ability to disperse these viruses to distant geographic locations. Coronaviruses (CoVs) repre- 
sent a family of zoonotic viruses with wide variety of animal hosts, including birds and 
humans. The infections caused by coronaviruses vary from mild to severe, depending on 
the viral species and the host. Since the coronaviruses exhibit extraordinary large RNA 
genome, also the rate of homologous recombination is high, which in turn contributes to 
the genetic diversity and interspecies host-switches of CoVs. The emergence of novel CoVs 
has been rich during the last decades, and wild birds seem to serve as reservoirs for a variety 
of CoV strains. We examined the CoVs circulating among wild birds in Finland. 

Materials and methods: Samples (cloacal swab, tracheal swab, oropharyngeal swab, or 
tissue) representing 61 bird species were collected during 2010-2013, and examined by RT- 
PCR targeting the RdRp gene for the presence of CoV RNA. 

Results: Altogether 51/939 (5.4%) of the examined birds were found positive by RT-PCR. 
Diverse gamma- and deltacoronavirus sequences were detected. 

Discussion: Gamma- and deltacoronaviruses circulate among wild birds in Finland. The 
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number of CoV-positive birds detected each year varies greatly. 


Introduction 


Coronaviruses (CoVs) comprise a family under the 
order Nidovirales (family Coronaviridae) and infect 
a wide variety of mammals and birds. The course of 
infection varies greatly from asymptomatic to severe 
disease, depending on the host and virus species in 
question. The genome of CoVs is one of the largest 
(25-32 kb) viral RNA-genomes [1]. Based on phy- 
logenetic analysis, the CoVs are divided into four 
different genera: Alpha-, Beta-, Gamma-, and 
Deltacoronavirus. The alpha- and betacoronaviruses 
are carried by mammals, whereas the gamma- and 
deltacoronaviruses mainly infect birds, with few 
exceptions [2,3]. The large genomes, infidelity of 
the RNA-dependent RNA polymerase, and high fre- 
quency of homologous RNA recombination are the 
main factors contributing to the high genetic diver- 
sity of CoVs [4-6]. 

The first CoV, Infectious bronchitis virus (IBV), 
was identified in 1937 [7]. IBV mainly infects chick- 
ens, but may infect other bird species as well. IBV is 
highly contagious and affects the respiratory tract, 
gut, kidney, and reproductive systems, causing sub- 
stantial economic losses in the poultry industry [8]. 
Despite the global distribution of IBV, poultry in 


Finland remained free of clinical cases until April 
2011 [9] after which outbreaks involving several 
CoV genotypes have occurred in Southern Finland. 
The first human CoVs were identified in 1960s [10- 
12]. The human CoVs cause generally mild to moder- 
ate upper respiratory tract infections [13-15]. In 2003, 
a novel highly pathogenic betacoronavirus emerged in 
China, causing severe disease characterized by acute 
respiratory distress and it became known as severe 
acute respiratory syndrome (SARS)-CoV. The emer- 
gence of SARS-CoV [16] inspired virologists to more 
explore the highly divergent group of coronaviruses 
and their hosts, leading to the identification of a 
rapidly growing number of CoV species particularly 
in bats [2]. More recently, another highly pathogenic 
betaCoV infecting humans, the Middle East 
Respiratory Syndrome (MERS) CoV emerged in 2012 
with a case fatality rate of over 40% [17,18]. 
Migratory birds have the ability to facilitate the 
dispersion of microorganisms with zoonotic poten- 
tial. Wild birds have been associated with the ecology 
and dispersal of at least West Nile virus, tick-borne 
encephalitis virus, influenza A virus (IAV) and 
Newcastle disease virus (NDV) [19-21]. Since the 
discovery of IBV in 1937, it remained the only 
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known Gammacoronavirus for over 50 years, but the 
number has increased dramatically during the last 
10 years [5]. Thereafter, representatives of the genera 
Gamma- and Deltacoronavirus have been isolated 
from both wild and domestic birds including species 
from the order Anseriformes, Pelecaniformes, 
Ciconiiformes, Galliformes, Columbiformes, and 
Charadriiformes [22-24]. In this report we provide a 
description of CoV species circulating in wild birds in 
Finland. Altogether 939 samples representing 61 dif- 
ferent bird species were collected during 2010-2013 
and examined for the presence of CoV RNA. 


Materials and methods 
Ethical statement 


All active and passive surveillance samples tested in 
this study, including samples provided by hunters, and 
injured or diseased birds were initially sent to the 
Finnish Food Safety Authority Evira (Helsinki, 
Finland) for avian influenza testing as part of a 
national surveillance program. All hunters had appro- 
priate permits, and no birds were killed for research 
purposes. The hunted birds were shot during the 
annual duck-hunting season, from non-endangered 
species, and thereby do not require ethical approval. 


Sampling 


The samples were collected as part of a national 
active and passive surveillance program, coordi- 
nated by the Finnish Food Safety Authority Evira 
and University of Helsinki. The sample panels 
included (1) hunted clinically healthy birds (active 
surveillance), (2) hunted clinically diseased birds 
(active surveillance), and (3) birds found dead (pas- 
sive surveillance) or (4) clinically diseased (passive 
surveillance) (Table $1). The sample material used 
for RNA detection was cloacal swab, tracheal swab, 
oropharyngeal swab, or tissue. The swab samples 
were collected using nylon swabs and stored and 
transported in Universal Transport Medium (both 
by Copan International). Prior to RNA isolation 
from the swab specimen (stored at —80°C), the 
samples were centrifuged to remove solid particles 
and supplemented with additional antibiotics 
(streptomycin-penicillin). The tissue samples were 
homogenized and separated by centrifugation and 
supplied with antibiotics (streptomycin-penicillin). 
The sample panel from 2010 consisted of 343 sam- 
ples, the 2011 panel of 171 samples, the 2012 panel 
of 287 samples, and the 2013 panel of 138 samples. 
The samples were also screened for influenza A 
virus RNA [25]. 


Detection of coronavirus RNA 


QIAamp Viral RNA Mini Kit was used to extract the 
RNA from the samples. The CoV screening was per- 
formed using degenerate primers [26] to amplify a 
179 nt region in the polymerase (RdRp) gene (Orf1b). 
All positive samples were further characterized by 
amplifying a longer (608-610 bp) region in the 
RdRp gene using previously described primers (for- 
ward primer 5'-TGGGWTGGGAYTAYCC 
WAARTGYGA-3' and reverse’ primer  5’- 
GCATWGTRTGYTGNGARCARAATTC-3’) [23]. 
All PCR reactions were performed on a BioRad 
PTC-100 Thermal cycler or Applied Biosystems 
Veriti Thermal cycler by a one-step RT-PCR protocol 
using QIAGEN One-Step RT-PCR Kit (Qiagen, 
Germany). The reaction volume for PCR was 25 ul 
containing 5 ul 5x buffer, 1 ul dNTP mix, 1 uM both 
primers, 0.2 ul RNAse inhibitor (Applied Biosystems, 
20 U/ul), 2 ul RNA, and 1 ul Qiagen One-Step RT- 
PCR Enzyme Mix with thermal conditions as follows: 
reverse transcription at 50°C for 30 min, initial PCR 
activation at 95°C for 15 min, followed by denatura- 
tion step at for 30 s, annealing at 50°C for 30 s, and 
extension at 72°C for 30 s. The PCR products were 
analyzed on 2% agarose gels and gel purified using 
QIAquick gel extraction kit (Qiagen). 


Sequencing and phylogenetic analysis 


The RT-PCR amplicons were sequenced in both direc- 
tions at the Sequencing Unit of Institute for Molecular 
Medicine Finland, (Helsinki), using BigDye v3.1 
chemistry and run by ABI3730xl DNA Analyzer. The 
sequences were assembled in BioEdit and aligned 
using ClustalW [27-29]. The reference sequences 
were obtained through BLAST search [30]. and from 
GenBank [31]. Phylogenetic analyses were based on 
sequence alignments of 464 nucleotides (position 
100-547 of the RdRp gene, by IBV (AY392086) 
count) within Orflab. Sequence data from 46 of our 
samples was included in the analysis and the rest were 
left out because of inadequate sequence quality. The 
phylogenetic tree was constructed from Maximum- 
likelihood phylogenies and bootstrap values calculated 
by 1000 replicates. Gaps were handled by partial dele- 
tion. The evolutionary model GTR + G + I was used 
after evaluation of the best-fit model according to 
Bayesian Information Criterion. The alignments, selec- 
tion of evolutionary model and phylogenetic analyses 
were performed within Mega5 [29]. 


Results 
Coronavirus RNA in wild birds 


We screened in total 939 samples from 61 bird spe- 
cies (Table S1) during a time period of four years 


Table 1. Bird samples positive for coronavirus RNA. 


Species Clinical status (no.) 

Anas platyrhynchos Hunted clinically healthy (27) and hunted clinically 
diseased (1) 

Anas crecca Hunted clinically healthy 

Larus argentatus Found dead 

Chroicocephalus Found dead 

ridibundus 

Larus fuscus Found dead 

Columba sp. Hunted clinically diseased (1) and found clinically 
diseased (1) 

Cygnus cygnus Found dead 

Clangula hyemalis Found dead 

Total 


(2010-2013) using a traditional, conserved RT-PCR 
targeting a 179 fragment of the RdRp gene (Orflab) 
of all coronavirus lineages. Altogether 51/939 (5.4%) 
of the examined birds were found to be CoV RNA 
positive (Tables 1 and 2), of which 27 were found 
healthy, and 24 dead or diseased. CoV RNA was 
detected in eight species (Anas platyrhynchos, Anas 
crecca, Clangula hyemalis, Cygnus cygnus, Larus 
argentatus, Chroicocephalus ridibundus, Larus fuscus, 
and Columpa sp.). CoV RNA was most abundant in 
the samples from 2010 (11%) and 2013 (7.2%), 
whereas during 2011 (0%) and 2012 (0.7%) only few 
of the tested samples were found positive. Influenza 
A virus (IAV) was also detected in four of the CoV 
positive samples (Table 1) [25,32]. The sampling 
locations are illustrated in Figure 1. 


Phylogenetic analyses of the detected avian 
coronaviruses 


The samples positive in the initial RT-PCR screen 
were consequently amplified for phylogenetic pur- 
poses by the method described by [23]. The number 
of reference strains for the alignments was unfortu- 
nately limited as GenBank deposited sequences vary 
in both location and length as a result of the wide 
array of approaches that have been employed for 
sequencing of coronaviruses. Table 2 shows the 
GenBank accession numbers of the viruses illustrated 
in the phylogenetic tree (Figure 2). 

The viruses detected from wild ducks shared up 
to 99% nucleotide identity with other known 
CoVs of the Gammacoronavirus genus. The phy- 
logenetically most closely related viruses are iso- 
lates from ducks in Russia and USA, China [33], 
Hong Kong [24] and South Korea [33,34]. The 
gammacoronaviruses that were detected in gulls 
in this study shared up to 99% sequence identity 
and clustered with previous isolates from gulls in 
Siberia and Alaska [23]. The viruses isolated from 
pigeons (Columba sp) (Fin6709 and Fin11782) and 
an isolate from a long-tailed duck (Clangula hye- 
malis) shared only up to 89% nucleotide identity 
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% (no. positive/negative) Gamma Delta  Co-infections (no.) 


21.7% (28/129) X HON2 (1) + H3N8 (2) 

16.4% (9/55) X H3N8 (1) 

9.6% (5/52) X 

4.2% (1/24) X 

8.3% (1/12) Xx 

3.6% (2/56) X 

3.6% (1/78) Xx 

1.3% (4/8) Xx 
51/412 


with the currently published strains, with the 
highest nucleotide identity with duck gammacor- 
onaviruses. Two isolates, AvCoV/Larus fuscus/ 
Finland/10059/2013. and AvCoV/Chroicocephalus 
ridibundus/Finland/10083/2013, share 83% and 
85% nucleotide identity with published avian and 
mammalian deltacoronaviruses, respectively. They 
clustered in the Deltacoronavirus clade in the ML 
trees, most closely to the White-eye CoV isolated 
in Hong Kong (JQ065044) (Figure 2). None of the 
viruses from wild birds in this study showed close 
relatedness to circulating IBV-strains. 


Discussion 


Several viruses, including zoonotic and economic- 
ally significant pathogens, are known to circulate 
among wild birds. CoVs are a large group of 
viruses infecting mammals and birds, including 
infectious bronchitis virus (IBV), which is a highly 
contagious coronavirus infecting chickens. A vari- 
ety of coronaviruses are known to circulate among 
wild birds. We set up this study to determine the 
occurrence of coronavirus infections in wild birds 
in Finland. Altogether 939 samples collected 
between 2010 and 2013 from wild birds were 
included in the study. The prevalence was fairly 
high (up to 11%) although with great annual varia- 
tion. The sample panels differed to some extent in 
terms of time and _ geographic distribution 
(Figure 2) and sampled species, which might partly 
explain the annual variation. For example, majority 
of the samples of 2010 were from duck species 
(Anas sp.) sampled for active surveillance purposes. 
In 2011-2013 on the contrary, most samples were 
derived through passive surveillance and included a 
smaller proportion of ducks. However, no firm 
conclusions can be drawn whether the observed 
difference in prevalence is due to some regular 
pattern, or due to the sampling methods (hetero- 
geneity of the sample material). 

Some of the CoV positive samples were also posi- 
tive for H3N8 and H9N2 IAV RNA. Similar double 
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Table 2. (Continued). 


GenBank accession number 


KX588631 


Location 


Date of sampling 


9 September 2013 
21 October 2013 


Clinical status* 


Species 


Strain 


Taxonomic Family of the host 


171 
171 
171 


3 


European Herring Gull (Larus argentatus) 
European Herring Gull (Larus argentatus) 
European Herring Gull (Larus argentatus) 


Pigeon (Columba sp.) 


Avian Coronavirus/Larus argentatus/Finland/10879/2013 
Avian Coronavirus/Larus argentatus/Finland/13125/2013 
Avian Coronavirus/Larus argentatus/Finland/12822/2012 
Avian Coronavirus/Columba sp./Finland/11782/2013 
Avian Coronavirus/Columba sp./Finland/6709/2012 


KX588633 


KX588634 


6 September 2012 
4 October 2013 
28 July 2013 


3 
2 


KX588632 


11 
11 


Columbidae 


KX588636 


Pigeon (Columba sp.) 


*Clinical status. 


#1 hunted and clinically healthy. 
#2 hunted clinically diseased. 


#3 found dead. 


#4 found clinically diseased. 
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Figure 1.A map illustrating the sampling locations. The 
sampling locations (sub-regional units, year 2013) are indi- 
cated in the Table 2, and are here illustrated by numbers on a 
map of Finland. 


infections have been reported by other researchers 
[35], and low pathogenic avian influenza viruses 
(LPAI) together with type-1 avian paramyxoviruses 
(APMV) are also commonly detected [36]. The co- 
infections observed here were in four apparently 
healthy ducks, three mallards and a teal, suggesting 
that infections by a number of low virulence viruses 
are well tolerated by migratory ducks. 

CoV infections in mallard ducks (Anas platyr- 
hynchos) have been documented in Sweden 
[26,35], and one of the studies reported the pre- 
valence of CoV infections to be 6.9% albeit with 
seasonal variation [35]. In a more recent study, 
Wille et al. documented the prevalence of CoVs 
among wild waterbirds in Sweden to be even 
18.7% [37]. A study from Norway reported the 
prevalence of CoV RNA among graylag goose 
(Anser anser) in Northern Europe to be as high 
as 38% in 2004, but they also documented great 
annual and geographical variation; in 2003 the 
prevalence was only 18% [38]. Also in England, 
wildfowl (Anseriformes) and waders 
(Charadriiformes) have been reported to carry 
CoVs [39]. A report from South Korea documen- 
ted a prevalence of CoV RNA to be 0.95% in 
wildfowl (Anseriformes), and they also observed 
annual variation in the detection of CoV RNA 
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a0) KP757890_Porcine_deltacoronavirus/GHN-AH-2004 
KJ768231_Swine_deltacoronavirus/USA/OhioCVM1/20 14 
= EF584908 Asian_leopard_cat_coronavirus/Guangxi/F230/2006 
FJ376622_AvCoV/Munia_coronavirus_HKU13-3514 
Ee) JQ065046_AvCoVMagpie-robin/HKU18 
77 FJ376619 AvCoV/Bulbul_coronavirus_HKU11-934 
FJ376621_AvCoV/Thrush_coronavirus_HKU12-600 
JQ065044_AvCoV/White-eye_coronavirus HKU16-6847 
3 @ Fin10083_AvCoV/Larus_ridibundus/Finland/100 83/2013 
el @ Fini0059_AvCoV/Larus_fuscus/Finland/10059/2013 
JX548304 AvCoV/Shorebird_coronavirus/USA/85/200 4 


100 @ Fin6709_AvCoV/Columba_sp/Finland/670 9/20 12 
@ Fin11782_AvCoV/Columba_sp/Finland/11782/20 13 
81 @ Fin14395_ AvCoV/Clangula_hyemalis/Finland/14395/20 13 
100) GU393337_IBV_serotype_lowa_97 
90 GU393336_IBV_serotype_Holte 
= GQ504722_|BV_strain_Georgia_1998 
AJ311317_IBV_strain_Beaudette_CK 
KPO 330 30_DdCoV/DK/Jiangsu/F64/20 14 
4, @ Fin8812_AvCoV/Anas_platyrhynchos/Finland/88 12/20 10 
@ Fin8802_AvCoV/Anas_platyrhynchos/Finland/8802/20 10 
3 KF437672_AvCoV/spot- billed_duck/isolate_12WB52 
KP033044_DdCoV/DK/Zhejiang/1407/2014 
KPO 33036_DdCoV/DK/Guangdong/H1421/2014 
@ Fin8790_AvCoV/Anas_platyrhynchos/Finland/8790/20 10 
97] @ Fin8782_AvCoV/Anas_platyrhynchos/Finland/87 82/20 10 
@ Fin8764_AvCoV/Anas_platyrhynchos/Finland/8764/20 10 
@ Fin8770_AvCoV/Anas_crecca/Finland/8770/2010 
JN788840_AvGoV/J1451/Anas_acuta/0 91230 
@ Fin13377_AvCoV/Anas_platyrhynchos/Finland/1 3377/20 10 
@ Fin8788_AvCoV/Anas_crecca/Finland/8788/20 10 
KF437675_AvCoV/spot-billed_duck/isolate_12WB62 
JN788785 AvCoV/J0554/Anas_clypeata/091127 
JN788780_AvCoV/Anas_crecca/J0109/091106 
oa] @ Fin8629_AvCoV/Anas_crecca/Finland/8629/20 10 
@ Fin8585_ AvCoV/Anas_platyrhynchos/Finland/8585/20 10 
@ Fin8589_AvCoV/Anas_crecca/Finland/8589/20 10 
JN788831_AvCoV/J1395/Anas_acuta/091230 
8 GU396668_AvCoV/Pintail_coronavirus_PBA-10 
GU396670_AvCoV/Pintail_coronavirus_PBA-16 
KF437666_AvCoViisolate_12WB180 
fw KJ741881_AvCoV/duck/strain_OH38 
@ Fin9209_AvCoV/Anas_crecca/Finland/920 9/2010 
@ Fin11130_AvCoV/Anas_platyrhynchos/Finland/1 1130/2013 
GU396673_AvCoV/Pintail_coronavirus_PBA-124 
5' GU396672_AvCoV/Pintail_coronavirus_PBA-37 
| KP033070_DdCoV/DK/Anhui/F82/20 14 
| KM454473_DdCoV/DK/GD/27/20 14 
sor GU396687_AvCoV/Rock_sandpiper_coronavirus CIR-65824 
GU396675_AvCoV/Western_sandpiper_coronavirus_KR-28 
6, @ Fin10877_AvCoV/Larus_argentatus/Finland/10877/2013 


70 


@ Fin10879 AvCoV/Larus_argentatus/Finland/10879/20 13 


GU396679_AvCoV/Black-headed_gull_coronavirus_CIR-66183 
@ Fin13125 AvCoV/Larus_argentatus/Finland/13125/2013 
on fy @ Fin9211_AvCoV/Larus_argentatus/Finland/921 1/2010 
@ Fin12822_AvCoV/Larus_argentatus/Finland/12822/20 12 
GU396682_AvCoV/Glaucous-winged_gull_CIR-66002 
GU396674_AvCoV/Glaucus-gull_coronavirus_PBA-173 


Duck CoV Finland 2010 (22 viruses} 


KM093892_AvCoV/Porphyrula_alleni/isolate_7M/2011 
@ Fin4983_AvCoV/Cygnus_cygnus/Finland/4983/20 13 
JN788879_AvCoV/Dendrocygna_javanica/isolate_KH0 8-0 852/0 80506 
97] GU39667_AvCoV/Brent_goose_coronavirus_KR-70 
GU396690_AvCoV/Snow_goose_coronavirus_WIR-159 
EU111742_Beluga_whale_coronavirus_SW1 


100 


-—————__1 
0.2 


Alpha CoV 


Beta CoV 


Figure 2. Maximum likelihood phylogeny of avian coronaviruses detected in Finland in wild birds. A 464 nucleotide long 
sequence region of the RdRp gene within Orflab was sequenced and aligned with reference sequences: the analysis includes 
alignments of 46 viruses detected in this study and a diverse set of reference strains obtained through BLAST search. Maximum- 
likelihood analysis using 1000 bootstrap replicates was used to infer tree topology. The evolutionary model GTR+G + | was used 
after evaluation of the best-fit model according to Bayesian Information Criterion. Bootstrap support values exceeding 70 are 
shown next to the nodes. Strains described in this study are marked by dots. Branches have been collapsed where our viruses 
share high sequence identity, and the branch wideness is relative to the number of viruses included. 


[34]. Based on the studies from South Korea [34] 
and Northern Europe [38], and the results we 
obtained, it seems that large annual variation in 
the circulating RNA seems to be characteristic to 
avian coronaviruses. A study in China reported 
Gammacoronavirus and Deltacoronavirus RNA in 
12% of the screened, apparently healthy, wild 
aquatic birds (Anseriformes, Pelecaniformes, and 
Ciconiiformes) [24], which is in concordance with 
our results. Another study, conducted in the 
Beringia area (encompassing areas of Alaska and 
Siberia), reported Gammacoronavirus RNA in 
6.4% of the examined birds (Anseriformes, 


Pelecaniformes, and  Charadriiformes) [23]. 
Interestingly, some of the sequences found in the 
Chinese study [24] were identical to those found 
from the Beringia area [23]. Also the Finnish Duck 
CoVs reported in this study cluster very closely 
with duck CoVs of the Gammacoronavirus genus 
from Siberia and China, which are connected by 
migratory routes. The close genetic relatedness of 
these strains (as well as those found in South 
Korea and England) indicates that migratory 
birds do have a specific role in dispersing CoV 
to distinct geographical locations. However, the 
sequence analyses might be somewhat biased due 


to the limited availability of reference strains. This 
and previous studies have mainly been focused on 
duck and shorebird species of the Northern hemi- 
sphere, hence the involvement of other avian spe- 
cies and a wider global distribution of avian CoVs 
remains unclear. 

Interestingly enough our results show that the genus 
Larus (order Charadriiformes) can serve as a host for 
both genera, the Deltacoronavirus (species: Larus fus- 
cus) and Gammacoronavirus (species: Larus argenatus). 
Also, the fact that we detected a Deltacoronavirus in 
Chroicocephalus ridibundus is interesting, since in 
another study [23] Gammacoronavirus RNA _ was 
found from the same bird species. It was previously 
reported, that the bird species from the order 
Anseriformes can be infected with both delta- and gam- 
macoronaviruses [24], and we now demonstrate that 
the same applies for the order Charadriiformes. The 
abundance of the detected CoV sequences originated 
from ducks supports the hypothesis that the Duck CoVs 
would persist as endemic in duck populations causing 
subclinical infections similarly to LPAI and APMV. 

To conclude, we show that there is a wide range of 
gammacoronaviruses and also some deltacoronaviruses 
circulating in wild birds in Finland. In addition, we 
show that the number of CoV-positive birds detected 
each year varies greatly. Due to the high recombination 
rate of CoVs, new species constantly arise that are able 
to adapt to new hosts and ecological niches [5]. Since 
the number of CoV sequences detected from birds was 
so abundant, it would be interesting to study what types 
of CoVs are circulating in small Finnish mammals, 
especially those in close contact with birds. 
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